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SUMMARY 

Expressions  are  presented  for  the  load  distribution  on  a  represent¬ 
ative  group  of  plan  forms  in  sideslip  at  supersonic  speeds.  These 
expressions  were  obtained  by  the  application  of  lifting— surface  theories 
based  on  the  linearized  equation  for  compressible  flow.  Sketches  of  the 
load  distributions  are  included. 


INTRODUCTION 

In  three  recent  reports  (references  1,  2,  and  3)  the  variations  of 
rolling  moment,  of  yawing  moment,  and  of  lift  and  pitching  moment  with 
sideslip  have  been  investigated  for  a  group  of  wing  plan  forms  for  super¬ 
sonic  speeds.  The  pressure  distributions  required  to  compute  these 
forces  and  moments  were  calculated  using  linearized  compressible— flow 
theory  for  thin  airfoils.  Since  the  reports  referred  to  were  concerned 
with  the  detailed  expressions  of  moments  and  forces  for  the  various  plan 
forms,  it  was  decided  that  the  reference  value  of  the  pressure  distribu¬ 
tions,  their  possible  utility  in  stress  analysis  and  design,  and  the 
desirability  of  including  some  pictorial  representations  justified  the 
treatment  of  these  distributions  as  the  subject  of  a  separate  report. 

By  virtue  of  the  many  approximations  involved  in  its  derivation, 
the  linearized  theory  applied  constitutes  one  of  the  most  simplified 
analytical  approaches  to  compressible— flow  problems.  Furthermore,  in 
addition  to  the  factors  approximated  in  the  linearization  of  the  poten¬ 
tial  theory,  the  analysis  employed  does  not  account  for  the  lack  of 
complete  rigidity  of  a  wing  nor  the  effects  of  viscosity  in  the  flow. 
These  are  two  important  factors  that  may  have  considerable  effect  on 
the  actual  distribution  of  the  pressure  on  a  wing.  Thus,  it  is  not 
expected  that  these  pressure  distributions  will  conform  precisely  to 
those  obtained  in  the  actual  physical  flow.  It  is  expected,  however, 
that  these  theoretical  solutions  will  be  good  first  approxin£ations  for 
the  plan  forms  and  conditions  considered  herein  and  they  should  provide 
satisfactory  indications  of  the  pressure-difference  contours  in  general 
if  not  in  detail. 
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The  investigation  covers  the  following  configurations:  (See  figs. 

1  and  2.)  (l)  Triangular  plan  forms  with  subsonic  leading  edges  or  with 

supersonic  leading  edges;  (2)  trapezoidal  plan  forms  with  all  possible 
combinations  of  raked— in,  raked-out,  subsonic,  or  supersonic  tips;  (3) 
rectangular  plan  forms;  and  (4)  two  swept-back  plan  forms  with  super¬ 
sonic  trailing  edges  developed  from  the  triangular  wings.  Illustrations 
are  included  in  order  to  provide  a  convenient  visual  correlation  between 
the  expressions  for  the  pressure  distributions  and  for  the  moments  and 
forces  that  were  calculated  from  them  (references  1,  2,  and  3) ■  The 
arrangement  of  the  appendixes  was  based  on  the  desire  to  present  a  syste¬ 
matic  and  convenient  compilation  of  expressions  and  illustrations  for  the 
load  distributions  for  the  various  plan  forms  considered. 
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SYMBOLS,  COEFFICIENTS,  AND  AXES 

aspect  ratio 

span  of  wing  measured  normal  to  plane  of  symmetry 
Mach  number  parameter 

ratio  of  tangent  of  right  tip  angle  to  tangent  of  Mach 
cone  angle 


*/  M3.2— 1^) 


\  tan  p  / 


,  tan  p 

chord  of  wing  in  plane  of  symmetry 


incomplete  elliptic  integral  of  the  second  kind  with 
/  n<P 

modulus  k  ( 

complete  elliptic  integral  of  the  second  kind  with  modulus  k 


- - -  \ 

J  1-k2  sin2  9  d0  J 


E  (£,k) 


incomplete  elliptic  integral  of  the  first  kind  with  modulus  k 

0nCP _ d0_ 

n  V  1— k2  sin2  9 


) 


complete  elliptic  integral  of  the  first  kind  with  modulus  k 


F  (§,k) 


over-all  longitudinal  length  of  swept-back  wing 

slope  of  right  wing  tip  measured  in  plane  of  wing  (positive 
for  raked-out  tip,  negative  for  raked— in  tip) 
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Mx 

n0,nx 

AP 

AP 

q 

q 


free— stream  Mach  number 


slopes  of  plan— form  edges  relative  to  wind  axes 
pressure  differential  across  wing  surface,  positive  upward 


loading  coefficient 


free— stream  dynamic  pressure 


S 

u 

V 

w 

x,y,z 

a 

0 

P 

p 
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area  of  wing 

perturbation  velocity  parallel  to  positive  x  axis 
free— stream  velocity 

perturbation  velocity  parallel  to  z  axis  (positive  upward) 
rectangular  coordinates  of  wind  axes  (fig.  3) 
angle  of  attack,  radians 

sideslip  angle  (positive  when  sideslipping  to  right),  degrees 

s) 

rectangular  coordinates  of  stability  axes  (fig.  3) 
rectangular  coordinates  of  body  axes  (fig.  3) 
air  density  in  the  free  stream 
perturbation  velocity  potential 


Subscripts 

A,B,...V  expressions  given  in  Appendix  B 

The  body  axes  are  generally  a  right-handed  system  of  three  orthog¬ 
onal  axes  as  shown  in  figure  3  with  the  longitudinal  axis  £'  lying  in 
the  plane  of  the  wing.  The  stability  axes  are,  in  effect,  the  body  axes 
rotated  about  the  lateral  axis  tj*  (through  -a)  until  the  longitudinal 
axis  is  in  the  horizontal  plane  containing  the  free— stream  vector j  a 
subsequent  rotation  about  the  vertical  axis  £  (through  (3 )  would  bring 
the  longitudinal  axis  in  line  with  the  free— stream  vector  and  the  axes 
would  now  be  coincident  with  the  wind  axes. 
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It  should  he  mentioned  that  the  orientation  of  these  axes  as  shown 
in  figure  3  is  convenient  for  the  calculations  required  to  determine 
the  pressure  distributions  and  the  resulting  forces  and  moments.  In  the 
application  of  these  results  to  the  calculation  of  the  motion  and  dynamic 
stability  of  an  airplane.,  however,  the  axes  are  usually  rotated  so  that 
the  positive  direction  of  the  longitudinal  axis  is  into  the  free  stream 
and  the  positive  direction  of  the  vertical  axis  is  downward,  that  is, 
toward  the  undercarriage  of  the  airplane . 


METHOD  OF  ANALYSIS 

The  development  of  the  expressions  for  the  pressure  (i.e.,  load) 
distributions  on  wing  plan  forms  in  sideslip  was  merely  an  application 
of  supersonic  wing  theory.  In  this  report,  only  the  part  of  the  theory 
relating  to  the  flat— plate  or  so-called  "additional"  loading,  which  is 
the  loading  resulting  from  a  change  in  angle  of  attack,  will  be  consid¬ 
ered.  The  loadings  due  to  camber  and  twist  are  usually  assessed  inde¬ 
pendently,  and  the  Siam  of  these  two  loadings  is  often  referred  to  as  the 
"basic"  load  distribution. 

By  linearization  of  the  partial  differential  equation  for  compres¬ 
sible  flow  it  is  possible  to  develop  a  simplified  lifting— surface  theory 
for  thin  airfoils.  The  linearization  is  made  possible  by  the  assumption 
that,  for  thin  airfoils,  the  perturbation  velocities  induced  by  the  air¬ 
foil  are  small  relative  to  the  free— stream  velocity.  If  the  free— stream 
velocity  vector  is  parallel  to  and  in  the  direction  of  the  positive  x 
axis  and  if  cp  denotes  the  perturbation  velocity  potential  for  isen— 
tropic  flow,  the  linearized  partial  differential  equation  for  steady- 
state  conditions  at  supersonic  velocities  is 


d2(P  _  d2cp  _  d2cp 
cix2  Sy2  dz2 


=  0 


where  Mi  is  the  Mach  number  of  the  free  stream.  There  have  been  a 
number  of  methods  developed  that  provide  means  of  fitting  solutions  of 
■this  equation  to  the  boundary  conditions  of  thin-airfoil  theory  (e.g., 
references  4  through  9) •  The  results  to  be  given  herein  were  determined 
through  the  general  use  of  source— sink  and  doublet  distributions  (refer¬ 
ences  4,  5,  6,  and  9).  In  particular,  the  method  of  reference  6  was 
applied  to  cases  where  a  subsonic  tip  occurs  in  conjunction  with  a 
supersonic  leading  edge  or  tip;  whereas  the  load  distributions  for  all 
other  edge  and  Mach  cone  arrangements  were  calculated  by  application  of 
the  methods  summarized  in  reference  9* 

The  first  step  in  the  analysis  is  the  establishment  of  the  boundary 
conditions.  For  thin  airfoils  the  boundary  conditions  are  usually 
restricted  to  the  z=0  plane.  Thus,  if  the  local  angles  of  attack  at 
various  spanwise  stations  of  the  wing  Eire  specified  and  it  is  assumed 
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that  the  wing  is  coincident  with  the  z=0  plane,  the  boundary  conditions 
are  set.  Next,  the  expressions  for  the  loading  —  and  the  angle  of 

attack  a  are  formulated  in  terms  of  parameters  that  can  be  related  to 
the  potential  solutions  of  the  differential  equation  and  to  the  boundary 
conditions.  For  linearized  theory  these  relationships  are 


where 


and 


—  =  —  (if  u  is  for  the  upper  surface) 
q  V 


u 


Jz=0 


where 

w  =  —  ~ 

_  z=0 

Thus  the  problem  is  reduced  to  determining  q)  in  such  a  manner  that 
—  ^  is  equal  to  the  specified  local  angle  of  attack  at  every  spanwise 
station  of  the  wing. 

The  general  problem  of  specifying  the  angle  of  attack  and  of  solv¬ 
ing  for  the  resulting  velocity  potential  is  one  that  usually  requires 
the  solution  of  an  integral  equation.  (See  reference  9°)  For  cases 
where  the  edges  of  the  plan  form  are  supersonic,  however,  the  lack  of 
interaction  between  the  upper  and  lower-  surfaces  of  the  wing  permits  the 
problem  to  be  solved  by  a  distribution  of  sources  in  accordance  with  the 
local  slopes  of  the  plan  form  and  a  straightforward  integration  of  their 
potentials.  The  triangular  and  trapezoidal  plan  forms  with  supersonic 
edges  were  treated  in  this  manner.  Likewise,  wherever  a  subsonic  edge 
is  in  conjunction  with  a  supersonic  leading  edge  or  tip  a  straightforward 
integration  can  be  employed.  For  this  case  reference  6  provides  a  method, 
based  on  the  consideration  of  the  upwash  between  the  subsonic  edge  and  the 
Mach  cone,  whereby  the  usual  operations  involved  in  the  solution  of  the 
integral  equation  are  eliminated.  In  general,  however,  it  is  necessary 
to  go  through  rather  involved  procedures  to  calculate  the  load  distribu¬ 
tion  when  the  camber,  twist,  and  angle  of  attack  of  the  plan  form  are 
specified.  These  procedures  are  discussed  in  reference  9  wherein,  for 
conical— flow  conditions,  a  loading  element  is  used  to  set  up  the  inte¬ 
gral  equation  and  then  the  usual  integral-equation  techniques* are 
employed  to  solve  it. 

It  should  be  mentioned  also  that,  when  the  trailing  edge  is  sub¬ 
sonic,  an  additional  stipulation  based  on  some  physical  concept  for  the 
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flow,  such  as  the  Kutta  condition,  which  is  applied  herein,  is  needed  in 
order  to  eliminate  all  but  one  of  an  infinite  number  of  potential  solu¬ 
tions  that  will  satisfy  the  boundary  conditions*. 


PRESENTATION  OF  RESULTS 

All  the  plan  forms  and  conditions  investigated  in  references  1,  2, 
and  3  were  made  up  using  the  five  following  combinations  of  straight 
edge  s : 

1.  A  supersonic  leading  edge  in  conjunction  with  a  subsonic 

leading  edge 

2.  A  supersonic  leading  edge  in  conjunction  with  a  subsonic 

trailing  edge 

3.  Two  subsonic  leading  edges 

k.  A  subsonic  leading  edge  in  conjunction  with  a  subsonic 
trailing  edge 

5.  Two  supersonic  leading  edges 

The  expressions  for  the  load  distributions  on  these  five  combina¬ 
tions  are  given  in  Appendix  A  in  terms  of  the  wind-axes  notation. 

In  order  to  provide  an  easy  correlation  between  the  load  distribu¬ 
tions  and  the  aerodynamic  characteristics  of  the  plan  forms  presented 
in  references  1,  2,  and  3,  Appendix  B  contains  the  expressions  for  the 
load  distributions  on  various  sectors  of  the  plan  forms  in  terms  of  the 
plan— form  parameters  and  the  body-axes  notation.  Since  the  plan  forms 
are  restricted  to  the  z=0  plane  for  the  purpose  of  analysis  by  the 
thin-airfoil  theory,  the  body-axes  notation,  in  a  sense,  refers  to  coor¬ 
dinates  on  the  projection  of  the  plan  form  onto  the  z=0  plane,  which 
corresponds  to  giving  the  coordinates  in  terms  of  the  stability-axes 
notation.  This  slight  ambiguity  between  the  body  axes  and  the  stability 
axes,  caused  by  the  assumptions  employed  in  thin-airfoil  theory,  should 
not  be  allowed  to  cause  any  doubt  about  the  direction  of  the  normal 
force.  This  force  acts  perpendicular  to  the  plate  and  in  a  direction 
parallel  to  the  £'  axis,  not  the  C,  axis. 

The  order  of  presentation  of  the  plan-form  sectors  in  Appendix  B 
of  this  report  is  a  duplicate  of  the  arrangements  of  the  Appendixes  B 
of  references  1,  2,  and  3-  A  sketch  of  the  load  distribution  is  pre¬ 
sented  for  each  sector  in  Appendix  B  in  order  to  provide  a  convenient 
visual  interpretation  of  the  load  distribution. 

Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 

Moffett  Field,  Calif.,  Oct.  10,  19^9. 
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APPENDIX  A 

SUMMARY  OF  EXPRESSIONS  FOR  LOAD  DISTRIBUTION  ON  WING 
ELEMENTS  IN  TERMS  OF  WIND-AXES  NOTATION 

Expressions  Apply  to  Crosshatched  Plan-Form  Areas 

1.  Supersonic  leading  edge  in  conjunction  with  a  subsonic  leading 

edge: 
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2.  Supersonic  leading  edge  in  conjunction  with  a  subsonic  trailing  edge: 


AP 

q 


(ni+nQ)  -  2ninc 


y 

-n0+  — 
u  x 


ni 


_  Z 
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where 


E  is  the  complete  elliptic  integral  of  the  second  kind  with  modulus 
V  1-G2 

l-B^no—  J  (l-B2ni2)  (l-£2n02) 


B(ni-no) 

4.  Subsonic  leading  edge  in  conjunction  with  a  subsonic  trailing  edge 


y/.i  \  L  =  i 

x  B  X  ni  .  X  F 

x  x  ~  “o 


AP  _  4a  /  x  no 

q  =  B  7  ni  _  I 

x 


when  n0  =  0 


when  ni  =  — 
B 


P  =  -_y  1-*/  1-B2ni2 

E  aT2 


P  =  i/C2Z 

Jt  V  1-Bn0 


when  nQ  >  0;  ni  <  — 


E  E(cp,k)  -  K  E  (<P,k) 
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k*  = 


k  = 

cp  = 


1+B2nin0-  J  ( 1-B2ni2 )  ( 1-B2n02 ) 

B(ni+no) 


Gi-Bno 

1-BGmo 


K  =  F( 


V  l-k'2 


E  =  E( 


sin-1 


J  Gi2-k'2 
Gik 


Two  supersonic  leading  edges: 


ro|a  ro|a 
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APPENDIX  B 

SUMMARY  OF  EXPRESSIONS  FOR  LOAD  DISTRIBUTION  IN 
TERMS  OF  BODY-AXES  NOTATION 

Expressions  apply  to  crosshatched  and  heavily 
shaded  plan— form  areas 

1.  Triangular  Wings: 

A. 


n* 


where 


E  is  the  complete  elliptic  integral  of  the  second  kind  with  modulus 
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when  m  <  tan  P  < 


1-Bm 


Pi  = 


/  [B(m+  tan  pj-Gi(l-m  tan  p)](l+m  tan  (3) 

[B(m-  tan  P)+Gi(l+m  tan  p)](l-m  tan  P)(l-Gi2) 

_ Gx+k* _ 

—  +  E  +  k *  {e  f  (cp,k)-K  E  (<p,k)] 


Gi2-k»s 


where 


(l-ffl2  tan2P)-B2(m2— tan2p) 

2B(l+m2)  tan  P 

//[(l+m  tan  P)2-B2(m-tan  P)g3[(l--m  tan  p)2-B2(m+tan  p)2] 

2B(l+m2)  tan  P 

k»  _  Gi(l+m  tan  P)+B(m-tan  p) 

(l+m  tan  p )  +GiB ( m-tan  p) 


2j.|2 


cp  =  sin-1 


Gx  -k 


K  =  F(-,k);  E  =  E(-,k) 
2  2 
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8a(m+tan  p) 


B2(m+tan  P)2-(l-ffl  tan  P)2 


[B(m+tan  P)-(l-ni  tan  P)  ](m-tan  0) 
[B(m-tan  P)+(l+m  tan  p)](m+tan  P) 
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E*.  tan  (3  =  —  -  —  (See  footnote  a.) 
B+m 


( 


4a^m2+2mB— l) 


£P\ _ 

^  “Ea*  it  a/b2  ( m 2  +2mB— 1 )  2-(  B  +2m-m2B ) 


n 

2  + 


sin-1 


B2(m2+2mB-l)[(mB-l)  +  (B+m)  ]-(B+2m-ffl2B)[  (B+m)-  (Bm~l)] 


B(m—  3_)[(B+m)2+(mB— l)2] 


+ 


^  /( m+2B-B2m)  —  lp  (B2+2mB-l) 


(B2+l)  (m+  lr) 


(S). . 


W  TT^+PmB-l  1 


e3*  J  B2  ( m2+2mB— 1 )  — ( B +2m-m2B ) 


aLeft  leading  edge  hits  Mach  cone  from  apex. 
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Swept— back  wing  components: 
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/  (m+  3— )(B+tan  p) 

_ i! _ 

(l-B  tan  P)—  (B+tan  P) 
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(f  ),.-■(?), 

*  (?),  -  (?), 

where  r| ■  in  should  be  replaced  by  +  (|  -  mcr) 

Tlr  in  ^  ^  should  be  replaced  by  ■q'  —  (—  —  mcr) 


B***. 


(f)  -(f)  *(?)  -(f) 

b2***  m  q  e 

where  ij*  in  should  be  replaced  by  tj'  +  —  mcr) 

M 

T)f  in  should  be  replaced  by  -q *  —  —  mcr) 


ro|  o' 
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where  ’  in 


\q  /N 


should  he  replaced  by  +  (5-  —  mcr) 


B2— tan2P 
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ZAP  \  _  (  AP\ 

4*  vi 

/  ap\  _  /  apN  +  /  ap\  _  /  ap\ 

\  q  /s2*  \  q  A  '  q  '  0  \  q  'S 


where  t\  ' 


should  he  replaced  by 


V 


should  he  replaced  hy  t^* 


b 

2 


S**. 


(f\.  if), 
(f)„Af)  if )it\ 

where  tj  •  in  f  ) 

V  ^  /N 

<•  '■  ("1 


should  he  replaced  hy  ■nf  +  ll 

i  2 

should  he  replaced  hy  rj '  —  — 


r 
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Figure  f  -The  triangular ,  trapezoidal ,  and  rectangular  p/an- 
form  types  investigated . 
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Brn^.1 


Figure  2 Swept-back  plan  forms  and  Mach  cone 
configurations  investigated. 


x}  y}  z  rectangular  coordinates  of  wind  axes 
£  y,  £  rectangular  coordinates  of  stability  axes 
r ectangu/ar  coordinates  of  body  axes 


Figure  3.  -  C  oordinate  axes  systems  used  in  analysis. 
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